DELAYED-MUTUAL COUPLING DYNAMICS OF LASERS: 
SCALING LAWS AND RESONANCES 
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Abstract. We consider a model for two lasers that are mutually coupled optoelectronically by 
modulating the pump of one laser with the intensity deviations of the other. Signal propagation 
time in the optoelectronic loop causes a significant delay leading to the onset of oscillatory output. 
Multiscale perturbation methods are used to describe the amplitude and period of oscillations as a 
function of the coupling strength and delay time. For weak coupling the oscillations have the laser's 
relaxation period, and the amplitude varies as the one-fourth power of the parameter deviations 
from the bifurcation point. For order-one coupling strength the period is determined as multiples 
of the delay time, and the amplitude varies with a square-root power law. Because we allow for 
independent control of the individual coupling constants, for certain parameter values there is an 
atypical amplitude-resonance phenomena. Finally, our theoretical results are consistent with recent 
experimental observations when the inclusion of a low-pass filter in the coupling loop is taken into 
account. 
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1. Introduction. In recent work, we presented experimental and simulation re- 
sults for two mutually coupled lasers ^ with time-delayed asymmetric coupling. The 
light emitted from one laser propagates through fiber-optic cables to a photodetector 
that generates an electronic signal proportional to the light-intensity deviations from 
steady state. The electronic signal may be attenuated or amplified before it modu- 
lates the pump current of the other laser. The propagation time of the signal in the 
optoelectronic path introduces a significant time delay, and the coupling strength in 
each direction can be controlled separately. Our work in ^ investigated how the time 
delay and asymmetric coupling led to oscillatory and pulsating laser output. In this 
paper, we present a more theoretical exploration of the dynamics that result from this 
coupling configuration. 

In [2] we investigated this same configuration (asymmetric pump coupling) but 
without including the effect of delay. In addition, the coupling constant from laser- 
2 to laser-1 was negative, while the coupling constant from laser-1 to laser-2 was 
positive. That is, we assumed that the electronic coupling signal from laser-2 was 
inverted before applying it to the pump of laser-1. For purely harmonic signals, 
having opposite-sign coupling constants is equivalent to time-delayed coupling when 
the delay is half of the period. We found that as the coupling constant from laser-2 to 
laser-1 is increased in magnitude: (i) There is a Hopf bifurcation to oscillatory output, 
(ii) For certain parameter values there exists a small-coupling resonance such that the 
amplitudes of both lasers are pulsating, (iii) For large coupling, laser-2 acts as a small- 
amplitude, nearly harmonic modulation to laser-1. Laser-1 exhibits period-doubling 
bifurcations to chaos and complex subharnionic resonances. 
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In this paper we include the effect of the delayed-couphng, which results in a 
coupled set of delay-differential equations (DDEs) . With delay there is again a Hopf 
bifurcation to oscillatory output. However, the delay allows for periodic oscillations 
not just at the laser's relaxation period as in [2], but at periods that are integer 
multiples of the delay; we refer to the former oscillatory solutions as "internal modes," 
and the latter as "external modes." We show that as the coupling strength is increased, 
the first instability is a Hopf bifurcation to the internal mode. Our nonlinear analysis 
then shows that the amplitude of the internal mode varies as the one-fourth root of 
the bifurcation parameter's deviation from the Hopf bifurcation point. 

However, in the experiments of ^ the output was not the internal mode but an 
external mode. The period of the oscillatory output was a multiple of the delay, and 
the amplitude of the oscillations varied as the square root of the bifurcation param- 
eter. The discrepancy is due to the fact that the optoelectronic coupling contains an 
intrinsic low-pass filter that attenuates the internal mode. We show that the filter 
selects the external mode with the greatest possible frequency passed by the low-pass 
filter. A nonlinear analysis then shows that the amplitude of the external modes does 
indeed vary with a square-root power law. 

For the coupled lasers without delay in P, we found that the parameters could be 
tuned to cause a resonance-type effect with respect to the coupling; more specifically, 
the coupling could be tuned to maximize the amplitude of the oscillatory output. We 
find that this resonance effect also occurs with the inclusion of the delay. However, 
there are some dramatic qualitative differences as the coupling parameters are varied. 
In the system with delay, the branch of periodic solutions in the bifurcation diagram 
can be smoothly folded to form parameter intervals of bi-stability. Then there is a 
critical value of the same coupling parameter beyond which the bi-stability disappears 
suddenly and non-smoothly. 

Our coupling scheme is an example of "incoherent coupling" jSj because it depends 
only on the laser's intensity and not on the complex electric field. This is because 
the intensity of one laser affects the other only indirectly; in our case, the intensity of 
one laser modulates the pump current of the other. In contrast, "coherent coupling" 
refers to when the optical field of one laser is directly injected into the cavity of the 
other laser. Analogously, a single laser with delayed and re-injected self-feedback 
would be called "coherent" feedback. Semiconductor lasers with delayed coherent 
feedback have been extensively studied because of the their widespread application 
in electronics and communication systems; there the laser's output signal may be 
reflected off external surfaces back into the laser. The Lang-Kobayashi DDEs ^ 
are the canonical model for semiconductor lasers with delayed self-feedback and have 
been used to investigate phenomena ranging from the onset of instabilities to the 
development of chaotic output referred to as "coherence collapse" (see for a review). 
Coupled sets of Lang-Kobayashi or related DDEs are most often used to described 
coherently-coupled semiconductor lasers with delay. Two recent studies jSJ |7] are in 
the same spirit as this paper as they track the number of and properties of oscillatory 
solutions that appear as the coupling strength is increased. 

The main simplification that results from considering incoherent coupling between 
the lasers is that neither the laser-frequency detuning nor the clcctric-field phase dif- 
ferences affect the dynamics of the coupled system Pieroux et al. [HI have shown 
that the electrooptical coupling we consider in this paper leads to an equivalent dy- 
namical model as incoherent coupling. There have recently been a number of other 
investigations of lasers with incoherent coupling. The main focus of many of these 
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works was on chaotic synchronization (see |1U| and included references). Two recent 
papers by Vicente et al. ^2 ^] consider a similar implementation of optoelectronic 
coupling to that we consider here, i.e., the intensity of one laser modulates the pump 
current of the other. Their investigation of the oscillatory solutions is mainly numer- 
ical and they discuss interesting phenomena such as amplitude quenching, frequency 
locking and routes to chaos. 

Delay-coupled relaxation and limit-cycle oscillators have been the subject of many 
investigations (see ^J^^EIEIEI^I and their included references). These systems 
are self-oscillatory and the amplitude is often fixed by the intrinsic properties of the 
oscillator. In many cases the system may be reduced to a phase equation, or a system 
of phase equations if the oscillators are coupled ^HI ■ 

In contrast, for laser systems the amplitude is an important dynamical variable. 
This is because the lasers we consider are weakly damped, nearly conservative sys- 
tems such that the steady state is the only asymptotically stable solution. Oscillations 
must be induced by external mechanisms such as modulation, injection or coupling 
|20|. Thus, the amplitude of the oscillations is highly dependent upon the external 
mechanism, in our case the coupling, rather than the individual laser. We should 
mention that for coupled limit-cycle oscillators there is the amplitude instability re- 
ferred to as "amplitude death," where for specific values of the coupling (and oscillator 
parameter values) the amplitude becomes zero. 

After presenting our model, we show that the linear-stability analysis predicts 
a critical value of the coupling constant (that depends upon the delay time) such 
that there is Hopf bifurcation to the internal mode, i.e., oscillations with the intrinsic 
relaxation period. We then use multiple-scale perturbation methods, modified to 
account for time delays [21], to analyze the long-time evolution of the internal mode. 
The results are a pair of complex Stuart-Landau DDEs for the oscillation amplitudes 
that include a slow-time delay term. We analyze the amplitude equations to determine 
the bifurcation properties of the internal mode. Of particular note is that we allow 
for independent control of the coupling constants; most other studies of coupled lasers 
and oscillators consider symmetric coupling, where the coupling is the same across 
all of the elements. The independent control is important because, as we will show, 
it allows for a singularity in the bifurcation equations that marks the "resonance" of 
large amplitude solutions. 

The Hopf bifurcation to the internal mode occurs for small coupling. However, 
as the coupling is increased, the external modes appear via Hopf bifurcations. We 
have extended our multiple-scale analysis to be able to describe the bifurcation of the 
external modes that occur for 0(1) coupling. The analytical challenge is that the 
delay terms remain a part of the leading-order problem. Describing the bifurcation of 
the external modes is important for comparing our results to those in experiments [22] , 
because the experimental system contains low-pass filters that attenuate the internal 
mode originating from the first Hopf bifurcation. 

An alternative to our multiple-scale method for deriving slow-time amplitude 
equations is to apply center-manifold theory with averaging I23L 1141 115|. However, 
the averaging method does not retain the delay in the slow time; that is, the ampli- 
tude equations are ordinary-differential equations, not DDEs. The distinction is not 
important for their or our investigation because we look for steady-state solutions of 
the amplitude equations. However, time-varying amplitudes may require considera- 
tion of the slow-time delay. 

Reddy et al. ^1 1^ considered delay-coupled Stuart-Landau DDEs similar to 
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the amplitude equations we derive here. However, the complex coefScients they use 
are appropriate for limit-cycle oscillators and not the lasers that we consider. They 
focus their work on the properties of synchronization and amplitude death. 

In the next section, we nondimensionalize the model for two-coupled lasers, with 
the result being the focus of the rest of the paper. We perform the linear-stability 
analysis in Sec. 01 Slow-time evolution equations for small and 0(1) coupling are 
derived in Sec. El and SecO respectively. We close with a discussion of the results. 

2. Class-B model. We consider two class-B lasers (201 |2S1 modeled by rate 
equations as 

^ = J- 1,2 

(2.1) '^^^,^^^A,{t)-(l + I,)D.j\, 

where Ij and Dj are the scaled intensity and population inversion, respectively, of 
each laser. Scaled time t is measured with respect to the cavity-decay constant k, 
t = (2k) tr, where measures real time, — 7j|/(2K) is a ratio of the inversion-decay 
rate, 7|| , to the cavity-decay rate, k, and measures the relative life-time of photons 
to excited electrons. A{t) is the dimensionless pump rate and corresponds to the 
energy input to the laser by an external source, e.g., another laser, an incoherent light 
source or an electronic current. The mass- action coupling term ID models "stimulated 
emission" ; a photon passing through the laser cavity stimulates an excited electron 
to drop to the lower energy level, resulting in one less excited electron and one more 
photon. Eqs. H2.1|l may be derived as a reduced model from semi-classical Maxwell- 
Bloch equations |2()l 1251 1^ : the latter use Maxwell's equations to describe the laser's 
electromagnetic fields coupled with the Bloch equations from quantum mechanics for 
the amplifying media. 

If the pump rate is a constant, Aj{t) — Ajo, then the laser relaxes to the steady 
state DjQ = 1, IjQ = Ajo — 1. To facilitate further analysis, we define new variables 
for the deviations from the non-zero steady state as 

(2.2) Ij ^ Ijo{l + yj), Dj = 1 + £jy/l~^Xj, s^ei^/hot. 

Our goal is to investigate the effects of coupling through the pump. In addition, we 
account for the effect of a delay when the signal from one laser takes a non-negligible 
time before affecting the other. Thus, the pump coupling is taken to be 

(2.3) Aj{t) = Ajo - IjnhVkit - Tk). 

Thus, we feed the intensity deviations ijk — {Ik — Iko)/Iko from the steady-state output 
of laser k to the pump of laser j; the strength of the coupling is controlled by 6k. The 
pump-coupling scheme allows for easy electronic control of the coupling signal. The 
negative coupling results from the configuration of the electronic coupling circuits 
in 1^. Finally, we assume that the decay constants of the two lasers are related by 
£2 = The new rescaled equations are 

-yi - exi{ai + byi) - S2y2{t - T2), 
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dyi 
dt 
dx\ 

~dt 



= (3x2(1 + y2), 

(2.4) ^ = /3[~y2 - ePx2{a2 + 6^2) - 5m[t - n)], 
where 

(2.5) fll = — 0,2 = , and6=vAo- 

V-'io -'20 

For notational convenience we have let s ^ t and dropped the subscript on ei(ei —t e). 
Rogister et ah j28| have considered a very similar model to ours; however, in their 
case the cross-coupling term is instantaneous, while the delay appears through self- 
feedback of each laser's own intensity. 

3. Linear stability of the steady state. 

3.1. Characteristic equation. In the new variables, the steady state is given 
by Xj = tjj = 0. The linear stability of the steady state is determined by studying the 
evolution of small perturbations, for which we obtain the characteristic equation 

(3.1) [A(A + eai) + 1][A(A + ea2/3') + - l3''5i62e-^^^'+^^^ = 0. 

The delay results in the exponential term exp(— A(ti -I-T2)). The transcendental form 
of the characteristic equation and, hence, the possibility of an infinite number of roots 
is typical for DDEs. 

If either 5j = 0, then there are only decaying oscillations, which indicates that a 
single uncoupled laser is a weakly damped oscillator; this is the general characteristic 
of "class-B" lasers For 5j ^ numerical simulations indicate there is a Hopf 

bifurcation as the 5j are increased. To identify the Hopf bifurcation point, we let 
\ = iuj and separate the characteristic equation into real and imaginary parts. After 
some algebra we obtain a single equation for the frequency and another equation for 
the value of A = ^1^2 at the bifurcation point. The equations are 

(3.2) = tan(wT^)[(l - lj'^){0^ - w^) _ aia2l3'^ lo"^] 

+ew(l-cj2)(ai +02/3^), 

(3.3) (3^^l = [{l-uj^){P^-uj^) 

-e^aia2f3''uj^]^ + [eLo{l - uj^)iai + a2/3')]^ 

where = ti + T2 is the "round-trip" delay time. The transcendental equation for 
ui can be solved numerically and its solution is then substituted into Eq. (|3.3I) to 
determine the value of the coupling at the bifurcation point A = Ah ■ 

To simplify further discussion we consider (3 = 1. This implies €2 = ^lylwjl^xii 
which is a specific relationship between the lasers' decay constants, 7|| .j and Hj, and 
the pump rate, AjQ\ if /3 = 1 and Aiq = A20 then the lasers are identical. The results 
are qualitatively the same for nearly identical lasers with j3 ~ \. To simplify notation, 
we define ci = ai + 02 and C2 = aia2 and have 

(3.4) = tan(wr^)[(l - u^f - e^C2U)'^] + ew(l - uj'^)ci, 

(3.5) Al = [(1 - u:^ f - ^C2U?f + [ec.(l - ^'')c^\\ 

For e ^ 1 the leading approximation to the frequency is given by = ±1 or a; = 
m7r/Ts, m an integer. We will refer to the former as the internal mode because 
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Fig. 3.1. Neutral stability curve (solid curve) for Hopf bifurcations described by 

Eq. The solid circle at w ^ I indicates the internal mode. The crosses (+) at uj mn/rs in- 

dicate the external modes. Because Ts = nn, then lu f» m/n for the external modes, (e = 0.01, ai = 
aa = 2,6 = l,/3 = i; 

this is the scaled laser-relaxation frequency; more specifically, an oscillatory solution 
of Eq. (|2.4|l with period w = 1 is called the internal mode. Similarly, oscillatory 
solutions with period lu = rmr/Ts, m an integer, are called external modes because 
their periods are determined by the delay. For any fixed value of the delay there is 
the internal mode and an infinite number of external modes. 

From Eq. we see that for almost all values of the delay, as the coupling is 

increased from zero, the bifurcation to the internal mode occurs at a lower value of 
the Ah than the external mode (this is true except when Ts = m:, n even, which 
we will discuss in later paragraphs). More specifically, for the internal mode with 
ix) K, 1 the Hopf bifurcation occurs when A// = O(e^) <C 1. For the external mode 
with Lo — 0(1) the value of the coupling at the Hopf bifurcation is A// = 0(1). This 
is illustrated in Fig. 13.11 where we have plotted Eq. ^6.b\ . We will determine more 
precise approximations for these bifurcation points below; however, we will first make 
additional qualitative observations. 

Let us fix the delay to be a multiple of tt, i.e., = n7r, n an integer. Then uj — ±1, 
corresponding to the internal mode, is an exact solution to the frequency equation 
Eq. (|3.4|) and A// = ±e^C2. The + is taken if n is odd and the — is taken if n is even 
(the sign association comes from the real and imaginary parts of the characteristic 
equation before uj and A were separated). In the present paper, we will consider 
5j > and A > 0; results for A < are qualitatively the same. Thus, we have the 
exact solution u — ±1, Ah — +e^C2 when Ts = mr, n odd. An analysis of Eq. (|3.5|l 
shows this occurs at a minimum in the curve A//(w) as illustrated in Fig. 13.11 

If Ts is not a multiple of tt, say mr < < (n + l)7r, the Hopf bifurcation to 
the internal mode still occurs before the bifurcation to any of the external modes. 
The frequency of the internal mode is w w 1 (instead of = 1), and A^f > e^C2 but 
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still O(e^). Thus, if we plot A// as a function of the delay t^, the minimum of the 
curves Ajj{ts) occur when = mr, n odd. This is illustrated in Fig. 13.31 and will be 
discussed further below. 

Now consider the external modes with lo ~ nnr/Ts. For negative couplings when 
Sj > (see Eq. (|2.3ll 'l and A > an analysis of the characteristic equation shows that 
m must be even. Thus, the external modes have frequencies u « 27r/rs, Att/ts-, etc. 
In Fig. I3.ll we illustrate three cases when Ts = mr, n = 1, 5, and 53. The external 
modes then have to — m/n, m even. 

Finally, we consider the case when Ts = mr, for n even. As discussed in the 
previous paragraph, for negative couplings the external modes with frequency uj = 
m,Ti/Tg require that m = 2, 4, 6 . . .. Thus, for each external mode there is a value of 
Ts = mr, n even, such that w ~ I similar to the internal mode. This is illustrated in the 
numerical data of Fig. 13.3b . for = 2tt, while in Fib. 13.3b we see that the bifurcations 
to these two modes occur at nearly the same value of the coupling. Referencing 
Fig. 13.11 the two modes would have frequencies on opposite sides of w = 1 such that 
Ah is the same for each. This is referred to as a "double Hopf ' |2] EH) and may lead 
to more complicated bifurcation scenarios that we will not pursue in the present paper. 
We also note that at the double Hopf bifurcation, the designation of "internal mode" 
versus "external mode" is ambiguous because both have a; « 1. Indeed, following the 
numerical data in Fig. 13.3b . as is varied through 27r, we see that n = 1 internal 
mode becomes the m — 2 external mode, while the m = 2 external mode becomes 
the n = 3 internal mode. 

3.2. Approximations to internal and external modes. As discussed above, 
for any delay when the coupling is increased, the first periodic solutions to appear 
correspond to the internal mode and will have frequency to ~ 1. For arbitrary delay 
we can find an approximation for the frequency by letting uj — 1 + euji + O(e^) in Eq. 
1)3. 4|l to find that uji satisfies 

(3.6) taTl{Ts)[AL0l - C2\ - 2ciLUl = 0, 

which is a simple quadratic for lui. The bifurcation point is then approximated as 

A% = e^[16u;t + i{cl - 2c^)ujI + c^], 

(3.7) =t'(\iJi+aX)(^iJi+al). 

We can obtain simpler results than Eqs. H3.6|l and H3.7|l if we require that Tg « riTr, 
n odd; that is, in Fig. 13. 21 we examine locally to the critical values w(rs) = 1 and the 
minimum of the curves A// (r^). If = wn + eTsi, then 

(3.8) ^ = i_,2^^^ + o(£2)^ A\ = t^cl + 0(e'). 

2ci 

For the external modes we let uj = ujq + euji + O(e^) with ujf) = mTr/rj, m even, 
and find that 



(3.9) UJ — ujQ 



UJqCi 



'{1-UJ^)Ts 

and the Hopf bifurcation point is approximated as 

(3.10) Al,^(l-^oY + ^— (l--o'r-^ + 
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Fig. 3.3. Same as Fia. I'i.sA but with a wider range in ui and A. The + indicate the numerical 
results. The thick curves are the analytical approximations of the internal modes (same as Ftp. I A'. 
The thin curves are the approximations of the external modes. 
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The external mode results are valid for Tg such that ljq ^ \; thus, Tg must be away 
from the double Hopf point when Tg = rmr. 

In Figs. l3.2l and l3.3l we plot the linear stability given by: (i) numerically evaluating 
Eqs. (|3.4|l and H3.5() (marked with +) and (ii) the asymptotic approximations from 
Eqs. (|3.6|I - H3.7|I and Eqs. H3.9|l - H3.10|l (solid curves) . Fig. l3.2l focuses on the bifurcation 
of the internal mode. For any fixed value of the delay time Tg, we increase A until 
the curve is crossed at A = Ah', the frequency is w « 1. In Fig. 13.31 we pan out so 
that both the internal and external modes are visible. From here we see that as A is 
increased, the bifurcation to the internal mode always occurs first. As A is increased 
further, we cross secondary curves indicating the bifurcations to the external modes. 
More specifically, suppose that t = n (n = 1). After the internal mode, the external 
modes with frequencies to = m/n = to/1, to = 2,4,... appear in sequence. For 
longer delay times, the order in which the external modes appear will be different; 
this can be seen in Fig. 13.31 for Tg w 11, where the bifurcation to the external mode 
with TO = 4 occurs before the mode with to = 2. Alternatively, in Fig. 13.1b when 
Tg — Stt, as A is increased the bifurcation to the periodic solutions occur in the order 
oj Ri 1, Ri 0.8, , ~ 1.2, , w 0.4, which corresponds to the internal mode, followed by 
the external modes uj — m/n, n — 5, to = 4, 6, 2, . . .. 

4. Hopf Bifurcation of the internal mode. 

4.1. Two-time scales. We use the method of multiple scales to analyze the 
oscillatory solutions that appear at the Hopf bifurcation points. For the uncoupled 
lasers, oscillations decay on an 0(e) time scale, which suggests that we introduce the 
slow time T = et; time derivatives become ^ = ^ + We analyze the nonlinear 
problem using perturbation expansions in powers of e^/^, e.g., Xj{t) = e^/^Xji{t,T) + 
exj2{t, T) + . . the relevant nonlinear terms and the slow-time derivative then balance 
at 0(e'^/^). Finally, for simplicity we set /3 = 1 indicating identical lasers; however, 
relaxing this assumption to nearly identical lasers (/3 = l+0{e)) does not qualitatively 
change the results. 

We now consider the effect of the two-time scale assumption on the delay term. 
With the additional slow time, the delay term becomes 

(4.1) y,it^T)^y{t-T,T-eT). 

If dj ^ 1, then consideration of the delay term is postponed until higher order and 
appears as part of the solvability condition for the slowly varying amplitude. However, 
if 5j — 0(1), then the leading-order problem will contain the delay term and analytical 
progress is much more difficult. 

If r ^ 0(l/e) (large delay), then er — 0(1) and the delay term must be retained 
in the slow argument. However, if t = 0(1), then the slow argument can be expanded 

as in] 

(4.2) y,{t - r) = v,(t -t,T)- ^T-^Vjit -t,T) + ... 

Now the delay of the slow argument is postponed to higher order. We shall see 
that this leads to simpler slow-evolution equations for (5^1 (see Eq. (|4.5|l '). and is 
necessary to make any progress at all for S = 0(1) (see SecEJ. 

It should be noted that care must be taken when using a series expansion of a 
delay term in a differential equation. The Taylor series may itself be justified but doing 
so can change the stability of the differential equation. A simple example is given in 
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pn| , while 132 provide more theoretical discussions concerning restrictions on the 
size of the delay. In our presentation we will check the validity of our approximations 
by comparing our analytical and numerical results. 

4.2. Bifurcation equation. From the linear-stability analysis, we know that 
the first bifurcation will be to the internal mode with w « 1 when A = 5i52 — O(e^) 
(we assume that Ts ^ nn, n even, and thus do not consider the case of the double-Hopf 
bifurcation). We present the case that both coupling constants are of the same order 
5j — edj. However, as long as A = O(e^) the bifurcation results are qualitatively the 
same. 

Proceeding with the multiple-scale analysis, we find that at the leading order, 
we obtain the solutions 

(4.3) yji{t,T) ^ Bj{T)e'* + c.c, Xji{t,T) = iBj{T)e'* + c.c, 

which exhibit oscillations with radial frequency 1 on the t time scale. To find the 
slow evolution of Bj{T), we must continue the analysis to 0{e^^'^). Then, to prevent 
the appearance of unbounded secular terms, we determine "solvability conditions" for 
the Bj{T). Due to the scalings of Sj and yj, the delay terms Sjyj{t — Tj,T — t^j) are 
0{e^^^) and contribute to the solvability condition. The final result is 

(4.4) ^ = -f i?,(T) - '-\B,iT)\'B,{T) + '-d^T - T,,,)e— ^ 

j,k — 1,2, j ^ fc, and where r^^k — e^fc. The effect of the delay in the slow time 
appears explicitly in the delay terms Bk{T — t^^). The delay in the fast time has 
resulted in the exponential terms e~''^'= . Eq. H4.4I) is valid for arbitrary values of the 
delay (we have not simplified the delay term Bk{T — r^^k))- However, according to 
Eq. (|4.2|l . if T = 0(1), then T^^k = 0(e) and delay in the slow time is postponed until 
higher order. This results in a simpler solvability condition where the delay in the 
slow time does not appear. 

Periodic solutions to the original laser equations, Eqs. H2.4|l . correspond to T- 
independent solutions to the solvability conditions. The conditions are the same for 
both Eqs. (|4.4() and (|4.5|l because the delay terms in Eq. H4.4() become constants. The 
full T-dependent solvability conditions, including delay, are required only to analyze 
the stability of the periodic solutions. 

To determine the amplitude and phase of the periodic solutions, we let Bj{T) = 
Rj{T)e^^^^'^\ define the phase difference ip = 62 — 9i and set the time derivatives to 
zero. We obtain 

— —aiRi — d2R2 sin{ip — T2), 
— — a2i?2 + diRi sin(V' + n), 

(4.6) - -liRl - Rl) + di^^ cos(V + Ti) - cos(V - T2). 

O 1X2 J^l 

We find that 

a2d2S2 sin(r2) - aidiSi sin(Ti) 



(4.7) tan(V') 



a2d2S2 cos(t2) + aidiSi cos(ti) 
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where 5*1 ~ and ^2 = i?2 ^^'^ amplitudes are found from the imphcit equations 



fljdjiS'j + a^dlSi + 2aia2did2 cos{ts)SiS2 — d\d\ sin^(Ts)S'i52 = 
did2 sin(rs)5'2(S'i - S2) + 3did2 cos(rs)(a2 - ai)5'iS'2 

(4.8) +i{aid\Sl - a2dlSl) = 0. 

Of particular note is that just as in the linear-stability analysis, the amplitudes of the 
periodic solutions depend on the sum of the delays = n + r2 . This corresponds to 
the effective round-trip time of information from laser-j returning to laser-j. 

In App. we present explicit solutions of Eq. 14.81 that specify Sj = as a 
function of the parameters and the delay Ts ■ The expressions are rather complicated 
but are easily evaluated numerically. To obtain simpler expressions, we consider 
Ts — {n + ^)7r, where ^ ^ 1 and n an odd integer. That is, we tune the delay time 
Tg to be near one of the minimums of the neutral stability curves in Fig. 13.21 The 
bifurcation equation is 

(4.9) i?2^3(ai+a2)^+e^2i + 0(C'), 
where 

(4.10) Di = ^ - 1, D2 = ^ - 1, 

and Z21 is specified in App. 1X1 

We compare the results of numerically computed bifurcation diagrams IBBl and our 
analytical results in Figs. 14. 114.31 In Figs. 14. II and 14. 31 = Stt and 537r, respectively, 
and hence f = in Eq. (|4.9p . The analytical bifurcation equation describes the 
numerical results quite well even for d2 far away from the bifurcation point. An 
expansion of Eq. (|4.9() near the bifurcation point shows that 

(4.11) i?2 ~ (A- Ah)'/*, A = 5i(52. 

We also note that Eq. (|4.9|l is identical to the bifurcation equation we derived in 
12] for two coupled lasers with opposite sign coupling. The coupled laser equations 
in are identical to Eq. (|2.4() except that in (i) there is no delay such that the 
coupling is instantaneous and (ii) we allowed the sign of one of the coupling constants 
to be positive. Near the bifurcation point, the small-amplitude periodic oscillations 
are nearly harmonic as indicated by Eq. H4.3|l . Thus, a positive coupling constant 
corresponds to an effective phase shift of half the period or, more generally, a phase 
shift of Ts = rnr (see App.IbI. Thus, we expect that when Tg is tuned to the minimum 
of the neutral stability curves in Fig. 13.21 the bifurcation equations of the two cases 
should be equivalent. 

In Fig. 14. 21 we show the case when C 7^ 0. Here we have tuned the delay Ts to be 
to the right of the minimum with Ts — {5 + 0.1)7r. The dashed curve is the result of 
numerically evaluating the full bifurcation equations, Eqs. (14.81) . The solid curve is the 
approximation Eq. (|4.9|l and it shows very good agreement with the numerical results. 
In general, we maintain good agreement as we tune away from the minimum up 
to Ts = {nzL 0.15)7r, n odd. The full bifurcation equation, Eq. Il4.8() . where we have 
not made any assumption on the delay, maintains good agreement out to at least 
Ts = (n±0.4)7r; beyond that however, it loses fidelity as we approach the double-Hopf 
bifurcations at Ts — (n± l)7r (t^ = n7r, n even). 
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Fig. 4.1. Internal mode, n = 5 with ti = T2 = n7r/2, and e = 0.01. Numerical data points 
indicated by *. The solid curve is the bifurcation equation determined by Eqs. lA.fXl - lA.fX for ^ = 0, 
or Eqs. lA.lV - lA.lX . Because Ts = ri + T2 = nn, ^ = so the bifurcation equations are equivalent 
(di = l,ai = aa = 2,6 = l,/3= Ij. 
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Fig. 4.2. Same as Fig. \4.1\ fn = 5) except that n = T2 = {n + 0.1)7r/2. Thus, Ts = (n + 0.1)7r, 
5 = O.Itt and the delay is tuned to the right of the minimum of the neutral- stability curve. The 
solid curves are the approximate bifurcation equations based on ^ <C 1 given by Eqs. IA.(H) - IA.^ . 
The dashed curves are the more general bifurcation equations given by Eqs. lA.lt - lA.St . All other 
parameters are the same as in Fig. \4.l\ 
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Fig. 4.3. Same as Fig. \4-l\ except that n = 53. 



Finally, we note that the first correction to the bifurcation equation Eq. (|4.9fl is 
linear in ^. Thus, the bifurcation equation is not symmetric about the minimums of 
the neutral-stability curves; that is, the amplitude will be larger on one side of the 
minimum than the other for a given perturbation of A above the minimun. 

4.3. Small-coupling resonance. The bifurcation equation Eq. (|4.9() is singular 
when its denominator is zero {D2 — 0), or 



(4.12) 



(A = (5i(52). If the parameters are such that the singular point is before the Hopf point. 
As < A/f, then the singularity can be ignored because the vanishing denominator 
will not occur when periodic solutions exist for A > A/f . However, if A5 > A^f, the 
bifurcation equation will exhibit the pole- type singularity. Near the singular point, 
the bifurcation equation predicts that the amplitude of the oscillations will become 
large corresponding to a resonance. The condition that A5 > A/f corresponds to 
S\ > e^a|. Thus, in laser-2 the increase in the population inversion provided by pump 
coupling is greater than the effective losses. In 2 , when the lasers were coupled 
without delay, we investigated this resonance in detail and the bifurcation equation 
matched the result of simulations very well. For the present case with delay coupling, 
however, we find that there are some dramatic differences. 

In Fig. 14.4b . we compare the analytical and numerical bifurcation branches. For 
di — 2.3, the bifurcation equation predicts a strong resonance peak. The numerical 
data for di — 2.3 (*) exhibits a small resonance peak and there is good agreement 
between the analytical and numerical curves both before and after the singularity. In 
Fig. 14.4b we increase di from 3.0 to 4.5 to follow the deformation of the bifurcation 
branch from an isolated resonance peak to a curve that forms a loop before continuing 
to higher values of c?i ; this folding of the bifurcation branch is new to the system with 
delay and was not observed in pj . We want to make clear that the bifurcation branch 
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is not intersecting itself. The apparent intersection results from projecting the full 
bifurcation curve of di versus (xi, yi, X2, 2/2)1 which does not intersect, onto the {di, X2) 
plane. 

Because the multiple-scale method is a local analysis and we assumed small- 
amplitude solutions, it is not unexpected that the analytical and numerical results 
match well near the bifurcation point. Similarly, in the vicinity of the resonance 
peak, the amplitude of the solutions becomes 0(1) and so it is not surprising that the 
analytical results fail to match the numerical results. That said, in our previous work 
with coupled lasers without delay ,2 , the analytically determined bifurcation equation 
matched the numerically computed resonance peak throughout the parameter range, 
i.e., even for 0(1) amplitudes; that was a bit surprising. We can surmise that in 
the present case it is the folding of the bifurcation branch that causes the mismatch. 
It may be that if we continued our analysis to higher order and derived additional 
nonlinear corrections to the bifurcation equation, we might better describe the folded 
branch. This would be a non-trivial and tedious calculation; but without doing so our 
supposition remains speculation. In j5j we used a different asymptotic method better 
suited to analyzing periodic solutions with 0(1) amplitudes with only partial success. 
Its application to the present problem with delay would be quite complicated and we 
have not made the attempt. 

In Fig. 14.4b we increase di to 5.0 and see that there is qualitative change in the 
bifurcation branch; there is no longer a loop, and the bifurcation branch continues 
to negative values of di . In Fig. 14.4b we see that the change in the nature of the 
bifurcation curve apppears quite abruptly at a critical value of di . Typically, one ex- 
pects smooth folding and unfolding of bifurcation branches as a parameter is changed 
|34j . In the present case, we have not been able to isolate an interval of di over which 
such an unfolding might occur. To our knowledge, such a discontinuous unfolding of 
a branch of solutions in a bifurcation diagram has not been previously described. 

Finally, if di is increased further the bifurcation branch folds back and continues 
for positive d2 ; there also appear to be secondary bifurcations but we have not explored 
these in any detail. We mention that as the delay Tg increases, or as e increases, the 
effect of the singularity is diminished. Conversely, with reduced values of e, the 
bifurcation branch may exhibit more complicated turns and folds. It is intuitively 
understandable that increasing the damping (e) will diminish any type of resonance 
phenomena. The role of the longer delay is unclear. 

5. Hopf Bifurcation of the external modes. In this section, we describe the 
external modes that emerge via Hopf bifurcations as the coupling is increased. We fix 
Si — 0(1) and use 62 as the bifurcation parameter with 62 = d2o + eS22 + ■ ■ ■■ We again 
have the slow time T = et and the perturbation expansion of x and y in powers of e^/^. 
The multiple-scale analysis is more complicated because the 0(1) couplings result in 
a leading-order problem that contains the delay terms. Due to the two time scales, 
the delay terms are of the form 6jyj(t — Tj,T — eTj). To make analytical progress we 
need to remove the slow delay from the leading-order problem. This requires that we 
assume er, <^ 1 such that by using Eq. 1)4. 2|l the slow delay is postponed to higher 
order. We then have a restriction on the size of the delay such that r — o(l/e). Thus, 
we find that our results fit well with the external modes corresponding to the case 
Tg = 5tt, Fig. 13.1b . but are inaccurate in the case = 537r, Fig. 13.1b . Finally, to 
simplify the presentation, we consider the case of equal delays ti = T2 — r, such that 
Ts = 2r; analysis of the general case with unequal delays can be carried out in the 
same way. 
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Fig. 4.4. Internal mode, n = 1 with Ts = nir. (a) The solid curve is the predicted bifurcation 
branch from Eqs. lA.fXl when di is tuned such that the singularity in the bifurcation equation occurs 
after the Hopf bifurcation. The * are the numerical bifurcation data for di = 2.3. (b) Numerical 
bifurcation branches as they deform between di = 3.0 and di = 5.2. (c) There is a critical value 
of di , where on either side of the critical value the bifurcation branch has a loop and proceeds to 
higher values of d\, or bends back without a loop to negative values ofdi. All other parameters are 
the same as in Fig. \4.1\ 



5.1. Leading order. The leading order 0(e^/^) problem is 
(5.1) = L ■ Xi{t,T) - D ■ X,it-T,T), 

where 
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and Xi{t,T) 



( xn{t,T) \ 

yiiit,T) 

X2l{t,T) 

V y2i{t,T) J 



We look for oscillatory solutions of the form Xi{t, T) — UiB{T) exp{iojt), where B{T) 
is a slowly varying scaler amplitude (to be determined from a solvability condition at 
0(e^/^)). Ui is a vector that is determined by substituting our ansatz into Eq. (|5.1|l 
to obtain 



(5.3) = J -Ui where J = 





-1 





-6206^ 


1 


— ILU 















-1 
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J 



For a nonzero solution Ui, we require det J — 0. This results in the same condition 
obtained from the leading-order linear-stability problem; specifically, 



(5.4) 

Thus, we have that 
(5.5) 



-i2ujT 



0. 



m = even, positive integer. 



For later reference it will be useful to note that 

e"'^'^'' = 1, and er"^^ ^ ±1 = v„ 



(5.6) 

where i>m — +1 if fn/2 is odd (e.g. m — 2,6, . . .) and i^,, 
4, 8, . . .) is even. Finally, we find that 



-1 if m/2 (e.g. m 



(5.7) 



1 

icuui 

V "1 J 



where ui — Vm\ 

V 020 



(5.8) 



5.2. Second order. At 0(e) the problem is 
d 



dt 



X2{t, T)=L- X2{t, T)-D- X2{t -T,T) + F2, 



where F2 = 



( \ 



V 2:21^21 / 



Because the homogeneous problem is the same as the 0(e^/^) problem we can, without 
loss of generality, set the homogeneous solution to 0. The inhomogeneous term F is 
proportional to exp(i2a;i) so that the solution is 



(5.9) 



X2{t,T) = B{TfU2e 



+ C.C. 



where U2 is specified in App. |21 
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5.3. Third order . At 0(e^/^) we find the solvability condition that determines 
the slow-evolution equation for B{T). The 0{e^^'^) problem is 



(5.10) 
where 

(5.11) 



-X^it, T) = L ■ X^{t,T) - D ■ X3 {t -T,T) + F3, 
ot 



I -aixii - 522V2i{t - T,T) + 52oT-§fy2i[t - T,T) - -§pxii ^ 
a;i2yii + a;iij/i2 - 



V 



-a2X2i + 5iT-^yii{t -T,T)- ^X2i 

a^22y21 + 2:2iy22 " wV^^ 



I 



The vector .Pj has terms proportional to exp(ict;i) and exp(i2ijjt) and the former will 
lead to solutions of the form (J/s -I- Vst) exp(iwt). The secular terms V-^t must be 
eliminated to prevent unbounded solutions for large which implies that a solvability 
condition must be imposed on ^3. The solvability condition is formulated as follows: 
We look for a solution to Eq. (|5.1Q(I of the form Xz — U exp{iujt) and at the same 
time identify the terms F in F3 proportional to exp(ja;t). We then obtain a algebraic 
system of equations for the vector U as 



(5.12) 
where 

(5.13) 



O^J-U 



F = 



I {-iujai - 522UlVm)B - {iuj - 520TUiV,n)BT \ 

iuj{u22-l)\B\'^B- Bt 

— iLUa2UlB — {iLU — 5lTVm)BT 

iujui{u2A - u\)\B\'^B — uiBt 



For U to have a non-zero solution, the Fredholm Alternative requires that - F = 0, 
where V is the solution to ■ V = (the superscript H refers to Hermitian). We 
find that = (iti, iujui, 1, iuj), and the resulting condition for the amplitude B{T) 
is 



(5.14) 



dB ^ , , 

— = [pi +iqi)B + (pr, 



iqn)B\B\' 



where and are given in App. [0 



5.4. Bifurcation equation. To analyze the solvability condition given by 
Eq. we let B(T) ^ R(T)e^'^^'^'> to obtain 



(5.15) 
(5.16) 



ipi+PnR^)R, 



dR 
df 

dd 2 

of =11+ InR ■ 



The equation for 6 determines the frequency correction as a function of the amplitude. 
The bifurcation equation is determined by considering steady-state solutions to the 
equation for R, and we find that 



(5.17) 



R' = -PL = - 



\6i 



r2Cij2|l — uj^\ 
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Fig. 5.2. Same as Fig. 15. ll excevt m = 6. All other parameter values are the same. 



The onset of oscillations occurs when R = Q and determines S22 at the Hopf bifurcation 
point; this result matches exactly that obtained in the linear-stability analysis. 

In Figs. 15. ll and 15.21 we compare the bifurcation equation Eq. H5.17|l to the nu- 
merically computed result. In each figure we have Ts = 57r = 0(1) (see Figs. I4.l1 and 
14.21 for analysis of the internal mode). In Fig. 15.11 the external mode corresponds to 
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Fig. 5.3. Upper bound on the value of uj as a function of di giving a sign change for parameter 
r2 in Eq. FT^ . 

m ^ A and its direction of bifurcation is subcritical. In Fig. 15.21 we have m — Q and 
the bifurcation is supercritical. In each case, there is good agreement between the 
numerical and analytical result local to the bifurcation point, where the multiple-scale 
analysis' validity is strongest. 

The direction of bifurcation (super- or subcritical) in Eq. H5.17|l is controlled by 
the sign of the constant ^2, which is given in Eq. I)C.5|I . Analysis of r2 shows that 
> in the interval -^2/5 < lu < uj^idi), where ujz (di ) is shown in Fig. 15.31 (the 
lower bound is a zero of the denominator of r2, while the upper bound is the sole real 
zero of the numerator). Thus, for external modes with frequencies within this interval, 
the bifurcation will be subcritical; this is the case for the external mode in Fig. 15.11 
For all other modes r2 < and the direction of the bifurcation will be supercritical. 

For reference, in Fig. 15.41 we show the bifurcation diagrams of the internal mode 
(cj « 1) and the external modes {lu = m/n, m = A and m = 6) for the case n = 5; that 
is, we have combined the bifurcation diagrams of Figs. ETll5.ll and l5.2l The subcritical 
bifurcation for the external mode m = A folds to provide an interval of hysteresis. 
The bifurcation branches are projections onto planes so that the intersections of the 
curves are not relevant. The periodic solutions corresponding to the external modes 
are unstable as they bifurcate from the unstable branch of steady-state solutions. 
The internal mode is stable until « 59, when a period-doubling sequence to chaos 
begins. 

Finally, we note that for the external modes 
(5.18) i?- [A- Ah]1/2. 

The amplitude of the external modes varies as the square root of the distance from 
the bifurcation point. In contrast, the internal mode varies as 1/4 the power of the 
deviation. 

6. Discussion. We have analyzed the output of two mutually coupled lasers, 
where the light intensity deviations from steady state modulated the pump of the other 
laser. The coupling strength in one direction is held fixed, while we examine the effect 
of increasing the coupling strength in the other direction. The signal-propagation time 
through the optical fiber and the electronic circuit causes a delay leading to a model 
that is a system of delay-differential equations. Linear-stability analysis finds that the 
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Fig. 5.4. Numerically computed bifurcation diagrams, both internal and external modes, for 
n = 5 and Ts = nn. The left-most branch corresponds to the internal mode shown in Fig. \4-l\ 
periodic solutions are stable until d2 ~ 59, when there is a period- doubling bifurcation. The middle 
and right-most branches are the continuation of the external modes shown in Fig. I5.il (m = 4j and 
Fig. \5.'iA (m = &), respectively; periodic solutions are unstable along these branches (e = 0.01, di = 
l,ai =a2 = 2,b=l,l3=l). 



steady-state solution becomes unstable at a Hopf bifurcation; we call the resulting 
periodic solution the internal mode because it oscillates at the laser's relaxation- 
oscillation frequency. As the coupling is increased, subsequent instabilities occur with 
frequencies determined by the round-trip delay time — Jtitt/ts, m = 2, 4, 6, . . .; we 
call these periodic solutions external modes. 

Using a multiple-scale analysis, we derive bifurcation equations for both the in- 
ternal and external modes. We find that the amplitude of the internal mode increases 
with the 1/4-root of the deviation from the bifurcation point and is supercritical; i.e., 
R ~ -|-(A — AhY^'^ ■ The amplitude of the external modes increase with the 1/2-root 
of the deviation and may be super- or subcritical; i.e., R ~ ±(A — Ah)^^^- Both the 
initial instability and the bifurcation results depend on the product of the coupling 
constants A = Si52. For the analysis of the internal mode, we assumed that both 
coupling constants were of the same relative size, dj = 0(e). However, other scalings 
satisfy the Hopf condition, e.g., 6% = 0(e^/^) and 62 = 0{e^^^). We have found that 
this does not change the qualitative properties of the bifurcation. 

We have focused our analysis on the onset of oscillatory instabilities and just 
beyond. We have not considered the stability of the internal and external modes as 
the coupling is increased further. However, numerical simulations indicated period- 
doubling bifurcations of individual modes as well as multimode behavior. 

With independent control of the individual coupling strengths, we have observed 
an atypical resonance phenomena. As in our previous work on two mutually coupled 
lasers without delay |2j , there is an interval of the coupling parameter S2 over which 
the amplitude of the oscillations becomes large; we have referred to this effect as 
a "resonance." The bifurcation equation, Eq. (|4.9|l . can describe the amplitude of 
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the oscillations for S2 above and below resonance. A singularity in the bifurcation 
equation indicates the existence of the resonance, but Eq. (|4.9|l does not describe the 
amplitude within the resonance well. The parametric form of the singularity indicates 
that physically, the coupling term provides an effective negative-damping that cancels 
with the lasers self-damping and hence provides a resonance effect. 

What is different from the results in T is that the bifurcation branch folds with a 
change in the coupling parameter Si; in Figs. l4. j| we see that the resonance peak forms 
a loop in the {di,X2) plane. We also find that instead of a smooth unfolding of the 
loop, there appears to be a critical value of the parameter di where the folded portion 
of the bifurcation branch abruptly disappears. To our knowledge, such a discontinuous 
unfolding of a branch of solutions in a bifurcation diagram has not been previously 
described. We do not understand the mechanism or manner by which the abrupt 
change occurs, but this will be the focus of future work. 

The application of multiple-scale perturbation techniques to DDEs is, to our 
knowledge, a relatively recent development |5T]- For our analysis of the internal 
modes, the delay terms appeared only as part of the solvability condition; thus, the 
multiple-scale analysis was relatively straightforward. In contrast, for the external 
modes that bifurcate when the coupling constants are 0(1), the delay terms are 
included in the leading-order problem. However, by looking for periodic solutions 
we are able to continue the analysis and formulate a solvability condition. Finally, as 
mentioned in the text, multiple-scale expansion of the delay term in Eq. H4.2|l amounts 
to a Taylor series expansion that, while it may be justified in and of itself, can lead to 
erroneous results for the DDE [30]. Thus, it is important to compare our analytical 
results to those from numerical simulations as an important check of the work. 

As discussed in the introduction, the method of averaging has also been used 
to analyze the weakly-nonlinear characteristics of delay problems. Averaging and 
the multiple-scale technique will lead to a similar slow-time evolution equation for 
the amplitudes. However, the multiple-scale technique accounts for the delay in the 
slow time where averaging does not. Because we focused on the existence of periodic 
solutions, the slow-time delay is removed and both averaging and multiple-scales give 
equivalent results. However, stability of the periodic solutions, or the investigation of 
more complicated phenomena such as quasiperiodicity, would require the slow-delay 
from the multiple-scale analysis. 

Finally, we finish with a discussion relating the results of our analysis to experi- 
mental results observed in .1 and 22 . The coupling circuit in the experimental system 
has two important characteristics. First, the signal is inverted, which results in the 
negative signs that appear in front of the coupling constants in Eq. (|2.4() . Second, the 
circuit acts as a low-pass filter that suppresses coupling of the relaxation-oscillations. 
We discuss the effect of each of these below. 

The signal inversion of the coupling circuit results in negative coupling constants 
in Eq. H2.4|l . However, both the linear-stability and the leading-order bifurcation 
results depend on the product of the coupling constants A = 6162- Thus, local to 
the Hopf bifurcation it does not matter if both 6j are negative or both are positive. 
More generally, our results depend only on whether A is positive or negative, not on 
the signs of the individual coupling constants. This is effectively a symmetry result, 
because for small amplitudes the oscillations are nearly harmonic. This means a sign 
change is merely a phase shift. However, we have observed in numerical simulations 
that when the amplitudes become larger such that the intensity is pulsating, the 
symmetry is lost such that positive coupling results in different system output from 
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negative coupling. 

As mentioned above, the optoelectronic coupling circuit in pp and [52] acts as a 
low-pass filter on the coupling signal. We do not account for this in our model of the 
system (recently, Illing and Gauthier 35 have analyzed a DDE where they explicitly 
account for the bandlimited response of their feedback system). However, because our 
results address the linear and nonlinear dynamics of both the internal and external 
modes, we can make a comparison between theory and experiment. 

The low-pass filter characteristic of the experimental coupling circuit attenuates 
the high-frequency relaxation-oscillations, corresponding to the internal mode. The 
result is that the experimental system oscillates at one of the low-frequency external 
modes. This is illustrated in Fig. 16.11 where we indicate the experimentally observed 
external mode with a (o) for different values of the delay. The indicate the 
theoretical value of Ah(lo) for a bifurcation to an external mode. Linear stability 
predicts that as the frequency increases the coupling strength, Ah, required for a 
Hopf bifurcation decreases (see also Fig. 13.11 for small u). However, there is a filter- 
cutoff frequency ujc such that the relaxation oscillations and higher-frequency external 
modes are attenuated and only the external modes with frequencies uJe < t^c are 
observed. 

For example, in ^ we observed oscillations with period equal to the round-trip 
delay time; specifically, w = 27r/T = Ztt/ts corresponding to the external mode with 
171 — 2. This indicates that the filter-cutoff' frequency, ujc, is such that 2Tr /ts < coc < 
A'k/ts. If LJc was greater than A'k/ts then we would have expected to observe the 
external mode with m ~ A because the value of the coupling at the Hopf-bifurcation 
point decreases as oje increases, i.e., it will bifurcate at a lower value of the coupling. 

More recent experiments [21] confirm that with longer delays, external modes 
with TO > 2 are exhibited; this is shown in the lower two plots in Fig. 16.11 That 
is, as the delay increases, the external modes with larger m will be below the cutoff 
frequency, oJe = tott/ts < oJc- And it is always the mode with largest frequency, but 
still below cutoff, that is exhibited because Ajy(w) is least for that frequency. Taking 
all five plots together, the experimental data suggest that the low-pass filter cutoff 
frequency is uic « 0.05 because all of the observed modes have frequency less than 
uJc- Finally, we add that it was observed in that the amplitude of the oscillations 
followed a square-root power law as a function of the coupling. This is exactly as 
predicted by the bifurcation equation Eq. H5.17|l for the external modes. 

Our comparison between the theory and the experiment does have limitations. We 
do not know the detailed filter characteristics of the optoelectronic coupling circuit. 
It is known that the filter profile is most certainly not a step function but is instead 
frequency dependent; it depends on the properties of the optical cable, the electronic 
coupling circuit and the frequency response of the laser to the modulated pump. In 
addition, for large values of the delay the frequency difference between the external 
modes becomes very small and many external modes are excited at nearly the same 
value of the coupling. For these cases the linear stability theory may be insufficient 
to identify the mode that is observed in the experiments. 

Acknowledgments. I.B.S acknowledges the support of the Office of Naval Re- 
search. 

Appendix A. Amplitudes of periodic solutions. 
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Fig. 6.1. For each value of delay, linear stability predicts a set of external modes with frequen- 
cies u>e = mn/Ts, m even, indicated by +. The frequencies of the external modes have been plotted 
on the linear- stability curve Af{ (oj) as in Fig. lil.ll External modes with higher frequency bifurcate 
first because Aij(cj) is decreasing. For each value of the delay, the (o) indicates the specific external 
mode observed in experiment. All of the observed modes are less than a cutoff frequency due to the 
low-pass filter characteristic of the coupling circuit. 



where 

(A. 2) N = 3Q [aidld2 sin^ Tg — a1di{ai + 02) cosTg] 

— Said^dj sin^ Tg + 3aidid2{3a2 + ai) sinrg costs 
-\-6ala2di{ai + 02) sinrs, 
(A. 3) D = Q sin [ afdi — did^ sin^ + 2 0102^2 cos r^] 

234 2 2 222 

-\-d1d2 sin Ts — 4 aia2did2 cos Ts sin — 4 aia2d2 sin Ts 
+2 a^a2(ii cos — aldid2 sin^ + 2 a\a\d2 

(A. 4) Q = d\d\ sin^ — 4 0102^1^2 cos — 4 a^Cg, 
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and 
(A.5) 



2aidi 



{did2 sin^Tg — 2 0201 cosTg — QsinTg) i?2- 



When the delay is tuned to be near one of the minimums of the neutral stability 
curves, i.e., = rmr + ^, ^ <JC 1, we have 



(A.6) 
(A.7) 



^20 + ^^21+0(^2^, 
^1 



Z2Q = 3(ai + 02)- 



Z21 
D2 



3(£'i + 1) 

2dlDl 
d\d2 
a\a2 

02(^2 



[Di{al + 2>alai) + (aa - oi)(d? - afj] 



1, 
1. 



^20 



■^(^20/0^ + ^21) +o(f), 



(A.8) a; = l-e( iz2o + ^a2/Di ) -^e 



\z2i + \a2{Dr + 1) 



+o{e). 



Appendix B. Delay and opposite-sign coupling. 

Consider a negative coupling constant such that 



(B.l) 



If 2/2 is harmonic with frequency oj = 1 and given by 2/2 = Acos(t), then 

dxi ^ . , , 
— - +d2Acos{t), 
dt 

—S2Acos{t — mT), n an odd integer, 
(B.2) -62Acos{t-Ts), 

where Ts = mr is the total "round-trip" or system-delay because the other coupling 
constant, ^i, is negative. 

Appendix C. External modes: coefficients. 

C.l. Second order: 0(e). 



(C.l) 

and 
(C.2) 



U2 = 



( iw(2u22"l) \ 

W22 

w{2u2i — u\) 
U24 



U22 



{l-^io^y-6,620' 



U2A 



(1 - 4a;2)2 _ j^j. 
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C.2. Third order: 0{e-'/'^). 

(C.3) p( = ri[(5i522T - 2a)^(ai + 02)], Qi = riifmUiOjKai + 02)-^ + 2S22]} 

O20 



(C.4) pn = rir2Ti/m — ^1, qn = rir22w^ , 

1 ,5i(.5i+<52o+4i^')-(l + u?)(l-2a;2)(l-4a;2) 



2((5i(52ot2 +4u;2)' (1 - 4^2)2 _ j^^. 
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